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This paper presents an improved k-e model for low -Reynolds number turbulence near a 
wall. The near-wall asymptotic behavior of the eddy viscosity and the pressure transport 
term in the turbulent kinetic energy equation is analyzed. Based on this analysis, a modi- 
fied eddy viscosity model, having correct near-wall behavior, is suggested, and a model for 
the pressure transport term in the fc-equation is proposed. In addition, a modeled dissipa- 
tion rate equation is reformulated. We use fully developed channel flows for model testing. 
The calculations using various k-e models are compared with direct numerical simulations. 
The results show that the present k-e model performs well in predicting the behavior of 
near- wall turbulence. Significant improvement over previous k-e models is obtained. 

1. Introduction 

The k-e model is one of the most widely utilized turbulence models for various turbulent 
flows of engineering interest. Patel et al.W extensively reviewed two-equation models 
which can be integrated down to the wall. One of their conclusions was that the damping 
functions used in turbulence models, especially the one for the eddy viscosity, need to be 
further modified in order to improve model performance. In fact, as we shall see later, 
many existing k-e models do not provide correct near-wall behavior of the eddy viscosity. 
In addition, an asymptotic analysis of near-wall behavior of turbulence shows that the 
pressure transport term in the turbulent kinetic energy equation is much larger than the 
turbulent transport term. This near-wall behavior is also observed in direct numerical 
simulation of fully developed channel flows (Mansour et aU 2 l, Kim et alJ 3 l). However, in 
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existing k-e models, this pressure transport term is either ignored or included in a turbulent 
transport model. 

In this paper, we will first analyze the near-wall asymptotic behavior of the eddy viscosity 
and the pressure transport term in the ^-equation, and then, in sections 2 and 3, propose 
models according to their near-wall behavior. In addition, in section 3, we reformulate a 
model equation for the dissipation rate following an argument similar to that of Lumley.W 
In section 4, we present the results of the calculations with various k-e models and compare 
them with direct numerical simulations. 

To analyze the near- wall asymptotic behavior of the eddy viscosity and other turbulent 
quantities, we expand the fluctuating velocities and pressure in Taylor series about the 
wall distance as follows: 

«i = hy + ciy 2 + d 3 y 3 + ... 

U 2 = c 2 y 2 + d 2 y 3 + ■■■ 

( 1 ) 

u 3 = b 3 y + c 3 y 2 + d 3 y 3 + .. . 
p = a p + b p y + c p y 2 + d p y 3 + . . . 

where uj, u 2 and u 3 are the velocity components in the direction of x, y and z. y is normal 
to the wall and x ,z are parallel to the wall. The coefficients a p ,bi,c 2 , ... are functions of 
x, z and t , where t is the time. Using the continuity and momentum equations, Mansour 
et alJ 2 l showed the following relations between the coefficients, 

2c 2 = ~{ b\ t i + 63,3) 

a P ,i = ( 2 ) 

Op , 3 = 2 uc 3 

where () ,- represents a derevative with respect to Xi. The eddy viscosity is usually defined 
as 

= MUiJ + Uj,i) - f (3) 
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where { ) stands for ensemble average and k = {u t u,)/2 is the turbulent kinetic energy. 
For plane shear flows, we can write from Eq. (3) 


— (uv) 

VT = dUfdy 


( 4 ) 


and using Eq. (1), we obtain the near-wall asymptotic behavior of the eddy viscosity: 


v t~q^ = -(bic 2 )y 3 + (-{M 2 + cic 2 ) + 2{&ic 2 )(ci))y 4 + 0(y 5 ) (5) 

That is, near the wall ut is 0(y 3 ), because dU/dy is usually 0(1). Any reasonable eddy 

viscosity model should have this near-wall behavior. We shall see later that many existing 
models do not have this near-wall behavior. For later use, let us examine also the near- 
wall asymptotic behavior of the turbulent kinetic energy k and its dissipation rate e = 
u{ UijUij ). Using Eq. (1), we obtain the following relations for the k and e: 

k = — * — y 2 +((6jC,) + (b 3 c 3 ))y 3 + 0(y 4 ) (6) 

l = (b\) + (bl) + 4((6j ci ) + (6 3 c,))y + 0(y 2 ) (7) 


In addition, the pressure transport term in the ^-equation, II = —^( Uip it ), becomes (using 
Eq.s (1) and (2)) 

n = -2K{Mi} + {M 3 »y + 0(y 2 ) (8) 

while the turbulent transport term in the fc-equation, — ( kui) i , can be estimated as 0(y 3 ). 
Therefore, the pressure transport term is much larger than the turbulent transport term 
near the wall. 


2. Eddy viscosity model 

In this section, we will propose a model for the eddy viscosity using its near-wall behavior 
described in the previous section. The eddy viscosity model can be in general written as 


» i)i 

ut = c u i 


( 9 ) 
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where u' and £' are the turbulent characteristic velocity and length scale, respectively. 
Depending on the way to specify velocity and length scales, the eddy viscosity model 
can be a mixing length model, a one-equation ( k ) model or a two-equation (e.g. k-e) 
model. For example, in plane shear flows, Prandtl’s mixing length model specifies the 
characteristic velocity with £'dU/dy. For near wall turbulence, the Van Driest mixing 
length model further damps the length scale to y[ 1 — exp(— y+ /A)\ where y + = u T y/u 
and u r is the friction velocity. For more advanced mixing length models, see Baldwin and 
Lomax^, and King^. One-equation ( k ) models use k 1 / 2 as the characteristic velocity, 
which is determined by the turbulent kinetic energy equation. In two-equation models, 
e.g. k-e models, the length scale is usually specified by k 3 / 2 /e, where e is determined by a 
dissipation rate equation. In this paper we will concentrate on two-equation models, which 
are usually written as 

k 2 

V T — C^f^— (10) 

where C M = 0.09, and is a damping function. The form of the damping function is quite 
critical to the prediction of the mean flow field. In fact, the mean velocity field mainly 
depends on the eddy viscosity model. Therefore it is important for an eddy viscosity model 
to have a proper near-wall behavior. We have examined the near-wall behavior of eddy 
viscosity models based on various k-e model equations. The results are listed in Table 1, 
which shows that some of the k-e models do not have the correct near-wall behavior of the 
eddy viscosity: u t = 0{y 3 ). 

The quantity k 3 ! 2 je is usually considered as a characteristic length scale (or the size) of 
the energy containing eddies, £'. One expects that near the wall, the size of these eddies 
should be order of the wall distance 0{y). However, Eq.s (6) and (7) show that k 3 / 2 je is 
0(y 3 ). Hence, k 3 ^ 2 je is not an appropriate quantity to represent the length scale of the 
large eddies near the wall. However, we can define a variable e as 


e = e — v- 


dk/dxi dk/dxi 
2k 


( 11 ) 
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which has a nice property: e approaches e away from the wall and is 0(y 2 ) near the wall, 
according to the Eqs. (6) and (7). Therefore, fc 3 / 2 /e is 0(y) near the wall, and is a proper 
quantity to characterize the length scale of the large eddies. With this length scale, the 
eddy viscosity should be written as 

*T = <V,y (12) 

Now in order for vt to have correct near-wall behavior, the damping function f ^ must be 
0{y) near the wall and approaches 1 away from the wall. The damping functions used in 
various k-e models are listed in Table 2. If we consider the presence of the wall as the main 
effect on the eddy viscosity, then we may assume f ^ is mainly a function of y + (defined as 
u T y/u, where u r the friction velocity). The form of can be determined quite accurately 
if we know ut, k and e, for example, from the direct numerical simulations. We may also 
optimize the following simple form by numerical experiments: 

fn = 1 - exp(-Giy + - a 2 y + 2 - a 3 y +3 - a 4 y +4 ) (13) 

The optimal values for channel flows are a\ = 6 x 10 -3 ,ct2 = 4 x 10 _4 ,a 3 = —2.5 x 
10 -6 ,a4 = 4 x 10 -9 . This form, Eq.(13), does provide the required near-wall behavior. It 
can be further optimized using the direct numerical simulation data. 

3. Modeled k-e equation 

To complete the eddy viscosity model, we need the modeled equations for the turbulent 
kinetic energy and its dissipation rate. This section will analyze the near-wall behavior of 
the fc-equation and propose a model for the pressure transport term with a proper near- 
wall behavior. The equation for the dissipation rate is also reformulated with a formal 
invariant analysis. 

Let us start with the equation for the turbulent kinetic energy, 

kt + Ujkj = D„ + T + U + P — e (14) 
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where D„, T and II represent the transport of the turbulent kinetic energy due to the 
viscosity, turbulent velocity and pressure, respectively. P and e are the rate of production 
and dissipation of the turbulent kinetic energy. The terms on the right hand side of Eq. 
(14) are defined as follows: 

D„ = ukjj 
T = -(kuj) tj 

1 ( 15 ) 

n = --(puj) tj 

P — (uiUj ) Ui t j 

Using Eq.s ( 1 ) and (2), we obtain the budget of the ^-equation near the wall, 


Du — u {{b\) + ( 63 )) + 6v({biCi) + ( b 3 c 3 ))y -f 0(y 2 ) 
T = 0(y*) 

n = -2i'((4,c 1 ) + {hc: t ))y + 0(y 2 ) 

P = 0(y 3 ) 

e = 1 '({6?) + (i>I))+M(6,Ci) + {b 3 c 3 ))y + 0(y 2 ) 


(16) 


This budget shows that the term II is much larger than the term T, and II cannot be 
neglected if we want the ^-equation be balanced at the level of 0{y). However, existing 
models do not consider this term or simply combine it with the term T and model them 
as 

-(*«>),; = {“*,>} j ( 17 ) 

In this paper, we propose a model for the pressure transport term n which has a similar 
form to that of the standard turbulent transport model, but with a coefficient to ensure 
its correct near -wall behavior, Eq. ( 8 ). The proposed model form of II is 


n = { 


Co 


Vt 


hi 1 - ex p(-y + )j °k 


k,j } J 


(18) 
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where Co is an adjustable model constant. Its optimal value for channel flows is 0.05. In 
addition, in some existing k-e models, it is assumed that e = 0 at the wall. In that case, in 
order to balance the term D v , a nonzero artificial term D must be added to the fc-equation. 
The form of D for various k-e models is listed in Table 3. Finally, the modeled fc-equation 
becomes 


dk dk_ 
dt +Uj dx j 


d ,, vr.dk 

dxj o k dxj 


] + n + vt { 


dUj 

dxj 


dUj , dU x 

dx i ' dxj 


-e + D 


(19) 


In the present model, D = 0, since e is nonzero at the wall. 
The exact dissipation rate equation is 


e,t + Ujcj = D e v +T c + + PMD 


(20) 


where D e v , T f and II 4 represent the diffusion rate of the dissipation rate due to the viscosity, 
turbulent velocity and pressure, respectively. PMD stands for the entire mechanism of 
the production and destruction of the dissipation rate e. The terms on the right hand side 
of the above equation are identified as follows: 


T ( = -v(ui t kUi >k Uj) t j 

2v . , 

n = - — {P,kUj,k),j 

PMD = 2v((ui'kUj,k) d* {V‘k,iV‘k,j)')Ui,j < ^v( y UjUi i k)lj i^kj 

- 2v(u,i' k u jik uij) - 2v 2 (u iik] u t ' kj ) 

The term II £ is usually neglected or combined with the term T f and modeled as 


( 21 ) 


T = 



To model PMD, we define $ by 


PMD = -v* 
k 


( 22 ) 
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At the level of the k-e model, we assume ^ is a function of u, 1 / 7 % k, e, e, Uij and Uij k . 
Since 'J' is an invariant, it must be a function of the invariants that can be constructed 
from these quantities. Therefore we can write 


* = uurUijkUijkX) 

e ee 


where Rt is the turbulent Reynolds number k 2 /ue. We expand ^ in a Taylor series about 
u T U ltJ U, <} /e and uu T Uij k U t j k k/ee, and take only the linear terms. We obtain 


'k = V’o + V’ x ~~ T ~~zr + $2Vv T Uij k Ui y j k —_ 

e ee 


(23) 


where the coefficients rfio, rfii and ipt are in general functions of Rt. Finally, the modeled 
dissipation rate equation becomes 


de TT de d .. ur.de. _ . e rr rr _ _ ee _ 

m + Vi dT j = a^ I( " + + c 'f'V'T v ‘J v ‘-> - + E 


(24) 


where, C\ and C2 are the model constants, and f\ and fi are functions of R t . The term 
E in the present model comes from the last term in Eq. (23): 


E — U{J k 


(25) 


where we have taken ^2 = — 1. The form of E and C \ , C 2 , f\ and fi for various k-e models 
are listed in Tables 3 and 4. 

4. Fully developed turbulent channel flow 

Flows which have self-similar solutions are particularly useful for accurate model testing, 
because their solutions are independent of initial conditions. Therefore, we do not need 
to accurately choose the initial conditions for the k and e. In this paper, we use a fully 
developed channel flow for model testing. This flow is the simplest wall bounded turbulent 
shear flow with a self-similar solution. However, the complex features of the turbulence, 
for example, the effect of the wall on shear turbulence, are remained. In addition, the 
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k-e equations for the channel flow are exactly one-dimensional steady problems, numerical 
calculations will be easy and accurate. Recently, the measurements^ 7 ] well verified the 
direct numerical simulation of a fully developed channel flow.] 3 ] These data are available 
for comparison with model predictions. 

Let h be the half width of the channel, u T the friction velocity and Re r the Reynolds num- 
ber defined as u T h/v. Let U, k,e, ut and y be the non-dimensional quantities, normalized 
by tt r , u*/h,v and h, respectively. The modeled equations for the channel flow become 


i-y 

dy 6r l + u t 


d ( 1 ^s u T-\dk, ,dU x 2 l 

dy{~Rt. 7 ^ + ( ~ + + UT ('dy} faT r 


- e = 0 


_d 

dy 


{ 



i 

Re r 


„ , ee ,<PU ,2 1 

- C 2 f 2J + My,) 


= 0 


(26) 

(27) 

(28) 


where 


vj = C^f^Re,. 

c 

(dk\2 

\ dv > 


k 2 


e = e — 


■ dy ■ 

2k Re T 


fn = equation(13) 

r , 0.4 Re T k 2 

1 7 = 1 -exp — ( — 

J 1.8 L V 6e 

Co 

C fA l - ex p(-y + )] 

The boundary conditions axe simple. At the wall, 


) 2 ] 


U = k = 0 

dk 
dy 


( dk \2 
V dv > 


€ = 


2k Re T 


and at the center of the channel, 


dk de 
dy dy 


(29) 


(30) 


( 31 ) 


The numerical solutions with various k-e models are obtained using Patankar and Spald- 
ing’s method.! 8 ] We use 65 nodes stretched in the half width of the channel in the same 
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way as done in the direct numerical simulations: 

Vi = 1 ~ cos ( 2(n -^1) ) i = 1 ’ 2 ’- 65 ( 32 ) 

The main results from different k-e models for Re r = 180 axe plotted in figures 1-6. All 
the calculations are compared with the direct numerical simulation data. Figure 1 shows 
that the model of Jones and Launder^ (JL) underpredicts the mean velocity as well as 
the peak value of the turbulent kinetic energy. In figure 2, Chien’s modelf 10 ! performs 
better than the JL model, but it overpredicts the mean velocity near the center of the 
channel as well as the turbulent kinetic energy. In these two models, c = 0 at the wall is 
used as the boundary condition, so the dissipation rate near the wall cannot be correctly 
predicted. Lam and Bremhorstf 1 ^ use a nonzero boundary condition for e and have made 
some improvement for the mean velocity and turbulent kinetic energy compared with the 
results of the JL model, see figure 3. However, the mean velocity is still underpredicted 
near the center of the channel, and the dissipation rate near the wall is not correct. In 
figure 4, the model of Nagano and Hishida^ 12 ! presents a very good prediction for the 
mean velocity and shear stress, while the peak value of k is underpredicted. Their main 
modification to the JL model is a change in the damping function f ^ and the form of E. 
A zero dissipation rate at the wall is used. Figure 5 presents the results of the present 
k-e model, which shows the improvement in the prediction of all quantities, including the 
dissipation rate. The eddy viscosity profiles for various k-e models are shown in figure 6. 
Overall, the present eddy viscosity model has better behavior than others. In figure 7, we 
show some results of present model for several high Reynolds number flows. 

From the model testing, we conclude that the present k -e model has made a significant im- 
provement over previous k-e models according to the comparison with the direct numerical 
simulations. We find that the improvement is mainly due to the modified eddy viscosity 
model and the model of the pressure transport term in the ^-equation. The proposed 
dissipation rate equation also shows a better near-wall behavior than the previous ones as 
shown in figures 1-5. 
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Table 1 Eddy viscosity and boundary condition for e in various k-e models 


Model 

Vj' 

BC: 6 

JL 

0(y 3 ) 

0 

Reynolds 

0(y 5 ) 

By* 

LB 

0(y 4 ) 

u alk 

V 8y* 

Chien 

0(y 3 ) 

0 

NH 

0(y 4 ) 

0 

Present 

0(y 3 ) 

u aik 

U dy> 


Table 2 Damping functions used in various k-e models 


Model 

h 

h 

fi 

JL 

eX P(l + R,/5o) 

1.0 

1 — .3exp(— R 2 ) 

Reynolds 

1 — exp(— .0198.??*) 

1.0 

[1 - .3exp(— i? 2 /9)]/ M 

LB 

[1 — exp(— .0165i?*)] 2 

x(l + %*) 

l + (f) 3 

1 -exp(-i? 2 ) 

Chien 

1 — exp( — .0ll5y + ) 

1.0 

1 - .22 exp(— i? 2 /36) 

NH 

[1 — exp(— y + /26.5)] 2 

1.0 

1 — .3 exp (— i? 2 ) 

Present 

Eq. (13) 

1.0 

1 — .22exp(— R 2 /36) 


Table 3 Model terms 

in various k- 

■e models 

Model 

n 

D 

E 

JL 

0 


2-r(0) 2 

Reynolds 

0 

0 

0 

LB 

0 

0 

0 

Chien 

0 

2i >K 
y 2 

~jr exp(-.5y + ) 

NH 

0 


^r(W„)(0) 2 

Present 

Eq. (18) 

0 

-r(0) 2 


Table 4 Model constants in various 

k-e models 

Model 

c „ c x 

c 2 

(Tk cr t 

JL 

.09 1.45 

2.0 

1.0 1.3 

Reynolds 

.084 1.0 

1.83 

1.69 1.3 

LB 

.09 1.44 

1.92 

1.0 1.3 

Chien 

.09 1.35 

1.8 

1.0 1.3 

NH 

.09 1.45 

1.9 

1.0 1.3 

Present 

.09 1.45 

2.0 

1.3 1.3 


R t = K 2 jve, R k = y/Ky/u, y + = u T y/v. 


13 



-1.0 -0.8 -0.6 -0.4 -02 0.U 0.0 4.0 8.0 12.0 16.0 20.0 24.0 















-0.8 -0.6 -0.4 -02 0.0 00 4.0 8.0 12.0 16.0 20.0 24.0 



FIGURE 3. - COMPARISON OF k-E MODEL Willi DIRECT NUMERICAL SIMULATION IN A FULLY DEVELOPED CHANNEL FLOW. LAM AND BREMHORST; 

O DNS. KIM F.l AL . 
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FIGURE l|. - COMPARISON OP k-E MODEL Wl III DIRECT NUMERICAL SIMULATION IN A FULLY DEVELOPED CHANNEL FLOW. NAGANO AND HISHIDA: 

O DNS. KIM ET AL. 
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FIGURE 6. - COMPARISON OF EIH1Y VISCOSITIES WITH DIRECT SIMULATION IN A FULLY DEVELOPED CHANNEL FLOW. 
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